Gaussification of quantum states of traveling light beams in atomic memory 
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We propose and investigate a protocol for Gaussification of quantum states of traveling light 
beams in an atomic quantum memory that couples to light via quantum non-demolition interaction. 
The protocol relies on a periodic switching between two different QND couplings and the total 
coupling strength scales only logarithmically with number of Gaussified light modes. The present 
scheme can be used to prepare entangled states of two distant atomic ensembles and to purify and 
Gaussify noisy non-Gaussian entangled states of light while simultaneously storing the purified state 
in atomic memories. 
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I. INTRODUCTION 



Quantum information processing with continuous vari- 
ables (CV QIP) exploits the information encoding into 
variables with continuous spectra such as quadrature op- 
erators of light modes or collective atomic spin of meso- 
scopic atomic ensembles 0, H|. This approach exhibits 
distinct advantages, such as deterministic preparation of 
entangled quantum states of light, highly efficient homo- 
dyne detection enabling deterministic quantum telepor- 
tation [js,Jj] and a suitable interface between light and 
atoms 0t3|- The latter relies on the collective enhance- 
ment of the atoms-light coupling by means of strong aux- 
iliary coherent light beam and a mesoscopic atomic en- 
semble containing large number of atoms. Under certain 
circumstances, the atoms-light coupling can be treated as 
a quantum non-demolition (QND) interaction [H, Q that 
enables deterministic entanglement of two distant atomic 
ensembles Q and implementation of a quantum memory 
for light @. 

A central role in CV QIP is played by Gaussian quan- 
tum states whose Wigner function has Gaussian form. 
On one hand, these states can be easily prepared ex- 
perimentally using coherent laser beams, passive linear 
optics and squeezers, and on the other hand they ad- 
mit efficient theoretical description in terms of covari- 
ance matrices and displacement vectors. Recall that co- 
variance matrix 7 associated with an M-mode state is 



defined as 7^ = ({An,-, Ark}}, where Arj 



(r,), 



r = (xi,px, . . . ,xm,Pm), and Xj and pj denote the am- 
plitude and phase quadrature operators of the jth mode, 
respectively. The quadrature operators satisfy the canon- 
ical commutation relations [xj,Pk] — iSjk- Interestingly, 
the Gaussian states turn out to be extremal in certain 
sense (i"o| . Among all states with a given covariance ma- 
trix, several important quantities such as distillable se- 
cret key rate or certain entanglement measures are min- 
imized by a Gaussian state. This observation plays a 
crucial role in proofs of the security of quantum key dis- 
tribution schemes with coherent states and homodyne 
detection An essential ingredient in the proof of 

extremality of Gaussian states is a specific symplectic 
transformation that produces balanced superpositions of 



quadratures of all M input modes (To| . In particular, for 
the output mode where all inputs appear with positive 
weight we have, 



1 M 



u 3 > 



3 = 1 
M 



(i) 



Here x, p represent the amplitude and phase quadrature 
of the output mode. If all M input modes are uncor- 
rected and prepared in the same state p with zero dis- 
placement, (xj) = (pj) = 0, then in the limit M — > 00 
the output state converges to a Gaussian state with the 
same covariance matrix as that of the input states [To| . 
Operation ([1]) can be implemented by mixing light beams 
on an array of beam splitters [HI, [l3| . The Gaussification 
procedure can be modified by performing measurements 
on the other output ports of the beam splitters and con- 
ditioning on the measurement outcomes. This latter ap- 
proach forms a core part of the entanglement distillation 
schemes for continuous- variable quantum states 

In this paper we propose and investigate a scheme for 
Gaussfication of states of traveling light beams mediated 
by their interaction with atomic quantum memory. The 
procedure is inspired by protocol for storage of states 
of light in atomic quantum memory [6( and by scheme 
for coherent-state information concentration in atomic 
memory (l9| . The scheme involves repeated switching 
between two different types of QND coupling which en- 
sures that the total QND coupling strength between light 
and atoms grows only logarithmically with the number 
of Gaussified light modes M. Another appealing fea- 
ture of the suggested scheme is that the light beams need 
not be perfectly synchronized and can arrive at different 
times [19| and the resulting Gaussified state is stored in a 
memory. This protocol can be used to prepare an atomic 
memory in an arbitrary Gaussian quantum state. It can 
also be used to purify and distill entangled states of light 
such that the purified entangled state is transferred from 
light onto atoms and stored in a pair of distant atomic 
memories. 
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FIG. 1: (Color online) Proposed scheme for Gaussification 
of quantum states of M propagating light modes in atomic 
quantum memory. Each light beam (Lj) interacts with an 
atomic ensemble (A) and is measured by homodyne detec- 
tor (HD). The atomic state is then displaced by an amount 
proportional to the measurement outcome. The QND interac- 
tion between light and atoms is periodically switched between 
Hamiltonians (|10[) and with the help of a fast polariza- 
tion controller (PC) and magnetic field pulses (B). For more 
details, see main text. 

The rest of the present paper is organized as follows. 
In Sec. II we will briefly review the quantum atoms-light 
interface and describe our protocol. In section III we will 
investigate how this protocol can be used to map an en- 
tangled Gaussian state of light beams onto two distant 
atomic memories. In Sec. IV we will then demonstrate 
that our scheme can be used for entanglement purifica- 
tion of phase-diffused two-mode squeezed states of light. 
Finally, Sec. V contains conclusions. 

II. GAUSSIFICATION IN ATOMIC MEMORY 

We consider configuration where light couples to the 
collective (pseudo) spin of an ensemble of atoms, J. Its 
cartesian components satisfy the SU(2) commutation re- 
lations, 

[Jy,J z ]=U x . (2) 

The atomic spins are initially all aligned along the x axis 
(e.g. by optical pumping) such that (J x ) = FN a, where 
Na is the total number of atoms and F denotes the total 
ground state angular momentum of a single atom. For 
large Na, J x can be treated as a classical quantity and the 
operator replaced by its mean value in the commutation 
relation ([2) . We can then define effective atomic quadra- 
ture operators xa = Jy/y/ (Jx) a.ndpA = Jzj y (Jx), that 
approximately satisfy the canonical commutation rela- 
tions, [xa,Pa] — i- Imagine a vertically polarized signal 
light beam together with an auxiliary strong coherent 
horizontally polarized light beam propagating through 
the atomic ensemble along the z axis, see Fig. 1. Under 
certain circumstances, the off-resonant atoms light cou- 
pling is governed by a quantum non-demolition (QND) 
interaction [aQ, 

H = Hkxlpa, (3) 

where n oc ^/NaNl, and Nl is the number of photons 
in the auxiliary horizontally polarized light beam. Note 



that alternatively, the two polarization modes can be re- 
placed by two spatial modes with atoms placed in one 
arm of a Mach-Zehnder interferometer [3] ■ In the Heisen- 
berg picture, the QND interaction transforms the atomic 
and light quadratures as follows, 



-out 

x A 


= X A 


n out 




Pa 


= PA, 






x° L ut 


= X L , 






pT 


= PL 



The QND coupling with k = 1 can be combined with 
measurement of the output light quadrature p° L ut and 
displacement of the atomic quadrature p A ut , p°A Ut ~^ 
p° A ut + gp L Ut , where g = i is the displacement gain. In 
this way we can partly emulate beam splitter operation 
by QND coupling [20( and we obtain 

X°S* =XA+X L) pT = l(PA+ PL)- (5) 

A naive M-fold repetition of this protocol with M light 
modes, fixed interaction strength k = 1 and variable feed- 
back gain g k = l/(k + 1), yields 

M 

x A Ut ~ X A H~ ^""^ x L,j, 
3=1 

This closely resembles the target transformation ([1]) but 
the atomic state is additionally squeezed by a factor of 
1/VM + l. This squeezing which grows with M would 
eventually make the state of the memory very fragile and 
would enhance any imperfections in the feedback dis- 
placement applied to the squeezed quadrature pa- An- 
other drawback of this approach is that the total squared 
coupling strength, 

M 

i=i 

which is proportional to the total number of photons in- 
teracting with the atomic ensemble, will increase linearly 
with M, Kl ot — M . However, it is desirable to keep K\ Qt 
as small as possible because certain decoherence effects 
in the atomic memory increase with K\ Qt . The QND 
coupling constant k can be expressed as k 2 — dr/, where 
d is resonant optical depth of the atomic ensemble and 
7] is the atomic depumping factor pH . |22T | . For a given 
experimental configuration the optical depth d is fixed 
and k can be increased only by increasing the number of 
photons Nl in the auxiliary optical beam, which however 
simultaneously increases also atomic depumping due to 
scattering of photons. It was shown in Ref. [21( that this 
decoherence effect can be modeled as a lossy Gaussian 
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channel with added thermal noise. Covariance matrix of 
the atomic state after QND coupling further transforms 



Ja = (1 - V)1A + 2?77v 



(8) 



where 7 vac = diag(l,l) is covariance matrix of vacuum 
noise. Clearly, the attenuation and thermal noise addi- 
tion described by Eq. © decreases efficiency of atoms- 
light coupling. To minimize this detrimental effect it is 
crucial to keep the total squared coupling K^ ot as small 
as possible. 

If we just reduce the coupling strength while keeping it 
fixed, ft = 1/ sJM , then with appropriate gains we obtain 



19], 



™OUt 

X A 



XA + 



1 M 



1 M 



(9) 



This procedure achieves a fixed total squared coupling 
strength K^ ot = 1 independent of M , but the output am- 
plitude quadrature x A ut suffers from an uncompensated 
noise represented by the term x A on the right hand side 
of Eq. © . This noise would preclude faithful Gaussifica- 
tion of light states even in the asymptotic limit M — > oo 
and would limit the amount of entanglement and purity 
of the state that can be created in the memory via Gaus- 
sification. The noise could be reduced by a presqueezing 
of the atomic quantum memory achieved by a QND mea- 
surement of the xa quadrature with coherent probe light. 
However, squeezing xa ~ > %a/V~M would require QND 
coupling ft 2 = M — 1 , so we recover linear scaling of K? ot 
with M. 

We now present an alternative protocol that avoids 
the unwanted squeezing, residual noise or presqueezing of 
atomic memory and achieves logarithmic scaling of K^ ot 
with M. The Gaussification is still based on repeated in- 
teraction of the atomic memory with many copies of the 
light state. However, we adjust the value of the coupling 
constant km in each step by controlling the intensity of 
the auxiliary strong coherent laser beam. We also repeat- 
edly switch between two different QND couplings, 



and 



H2N-I = foK2N-lXL,2N-lPA, 



H 2 N = -frn 2 NPL,2NXA- 



(10) 



(11) 



Switching between the two Hamiltonians (fit)]) and (jTTJ) 
can be performed for instance by rotating the collective 
atomic spin of the atomic ensemble by a magnetic pulse 
and by rotating the polarization state of the light beam 
by a fast electrooptical modulator. The protocol is de- 
signed such that after M steps the atomic quadrature 



operators are given by 



XAM 



Pa,m 



Cm I xa 



M 

E 



XL, 



CmVM 



12) 



where Co is a free parameter and C m depends on Co and 
M, as we shall see below. Note that Cm specifies the 
overall squeezing of the atomic state after M steps of the 
Gaussification protocol while Co represents an additional 
squeezing of the initial atomic state. 

When analyzing the protocol we must distinguish odd 
and even steps. Consider first an odd step, j — 2N + 
1. The atomic and light quadratures are transformed 
according to the interaction Hamiltonian (jTUJ) , 



™out 

x A 

pT 

XT* 
PT' 



XA,2N + K2N+lXL,2N+l, 

PA.2N, 

%L,2N+1, 

PL,2N+1 — K 2 N+lPA,2N- 



(13) 



The output light quadrature p° L ut is measured and the 
atomic quadrature p A ut is displaced by an amount of 
92N+\PL Ut i where 32^+1 denotes the feedback gain. After 
2N + 1 steps the atomic quadratures thus read 



%A,2N+1 — XA,2N + K2N+l%L,2N+li 
PA,2N+1 = (1 — 92N+lK*2N+l)PA,2N 



92N+lPL,2N+l- 

(14) 



In order to preserve the structure of the transformation 
(fT2l) we must set 



K2N+1 



92N+1 



2iV 



yW+1' 



V2N+ 1 
2{N+l)C 2N [ 



which yields 



C 



2JV+1 — 



2N + 2 



C 



2N- 



(15) 



(16) 



2N+ 1 

We next proceed with an even step of the protocol, 
j = 2N + 2. In the Heisenberg picture, the coupling pT|) 
gives rise to the transformations, 



pT 



pT 



XA,"2N+l, 

PA,2N+1 + K-2N+2PL,2N+2, 
%L : 2N+2 — K>2N+2%A,2N+1> 
PL,2N+2- 



(17) 



We measure the output light quadrature x° L ut and dis- 
place the atomic quadrature x A ut by g2N+2X L ut ■ If we 
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choose 



K2N+2 



92N+2 



C 2N+lv W+2' 
V2N + 2 



2N + 3 



a 



2N + 1, 



(18) 



then we preserve the structure of the target transforma- 
tion (JT5J) with 



C2N+2 



2N 



2N + 3 



C2JV+1- 



(19) 



By combining Eqs. ([TBI) and (fT9")l we obtain a recurrence 
formula for C2N, 



C2N+2 
which yields, 



27V + 2 



v/(2iV + 3)(27V+l) 



C 



2AT, 



c v/2jvTT {2 m) c 



(20) 



(21) 



The ratio Cm /Co is plotted in Fig. 2. We can see that 
this ratio lies in a narrow interval, C < C M < V2C 
and it quickly approaches a fixed asymptotic value, 



li ^ M — /ZE 

M-Kx> C V 2 ' 



(22) 



This indicates that the squeezing of the atomic state 
re main s bounded and is saturated at a v alue of C M = 
y^n/2Co- In particular, if we set Co = \f2pK then the 
transformation (|12[) becomes in the limit of infinite M 
exactly equivalent to the target operation (TTJ). For any 
finite M we can choose Co such that Cm — 1 will hold. 
In this case, the initial atomic quadratures will appear 
in the superposition rescaled by the squeezing factor Co- 
The optimal choice of Co will generally depend on the 
intended application of the protocol and on M. 
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FIG. 2: The ratio Cm /Co is plotted as a function of M (blue 
filled circles). The horizontal dashed line indicates the asymp- 
totic value of the ratio which is equal to ~ 1.253. 




FIG. 3: Dependence of the total coupling strength K tot on 
the number of steps M of the protocol is plotted for Co = 1 
(blue empty circles). The dashed curve represents the best fit 
by a logarithmic function. 



The asymptotically constant value of Cm/Cq implies 
that the QND coupling strength scales as km oc Cq/VM. 
The total squared coupling strength K^ ot defined by Eq. 
(O thus increases only logarithmically with M. This is 
graphically illustrated in Fig. 3 where the dependence of 
K^ ot on M is plotted for Co = 1 together with the best 
logarithmic fit K? ot « 0.05 + 1.12 ln(M + 1). 



III. MAPPING OF GAUSSIAN STATES OF 
LIGHT ONTO ATOMS 

If the Gaussification protocol is applied to light beams 
prepared in identical independent Gaussian states with 
covariance matrix 7^ and zero displacement, then the 
protocol maps the Gaussian state of light into atomic 
memory. Assuming that the atomic memory is initially 
prepared in some Gaussian state with covariance matrix 
j A and zero displacement, the covariance matrix of a 
Gaussian state in atomic quantum memory after M steps 
of the protocol reads 



(M) 

Ta 



1 



M + l 



SaIaSj 



AI 



M + l 



where 



S A 




Sl = 



Cm 




(23) 



(24) 



represent the effective squeezing of the atomic and light 
contributions. In the limit of large M and by choosing 
Co = \[2pK such that limM-s-oo Cm = 1 the many copies 
of a Gaussian state of light are mapped onto a single copy 
stored in the memory, 



M 



lim 7 ^ M) = 



1L- 



(25) 



This procedure can be used to establish entanglement 
between two distant atomic quantum memories by map- 
ping parts of an entangled state of light into them. The 
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Gaussification procedure is applied locally to each atomic 
memory and sequence of light modes. We shall assume 
that both atomic quantum memories are initially in a 
thermal state with mean number of quanta n and that the 
light modes are prepared in a pure two-mode squeezed 
vacuum state I^tms) with squeezing constant r and co- 
variance matrix 



1L 



( cosh(2r) 


sinh(2r) 

V o 





cosh(2r) 


— sinh(2r) 



sinh(2r) 


cosh(2r) 





The entanglement between two quantum memories can 
be conveniently quantified by a logarithmic negativity 
-E/v, which is an analytically computable measure of en- 
tanglement for two- mode Gaussian states 23]. We have 



E N = max(0, - log 2 (/x)), 



(27) 



where fj, is the lower symplectic eigenvalue of a two-mode 
covariance matrix corresponding to a partially trans- 
posed state of the two atomic memories, 

2 _ M 2 e- ir + M(2n + 1) (C 2 + Co" 2 ) + (2n + l) 2 
At ~ (M + l) 2 ' 

(28) 

The entanglement created in the atomic memory is max- 
imized when the symplectic eigenvalue /i is minimized. 
This occurs for Co = 1 and we have 



fJ-opt 



Me 



-2r 



2n 



M + l 



(29) 



The dependence of on n is plotted in Fig. 4 for 
various numbers of steps M of the protocol. We can 




FIG. 4: Entanglement of two atomic ensembles created by lo- 
cal mappings of parts of pure two-mode squeezed vacuum onto 
the atomic memories. The logarithmic negativity En of the 
two-mode atomic state is plotted as a function of the initial 
mean number of thermal quanta n in the memory. This de- 
pendence is shown for several number of Gaussification steps 
M — {1, 2, 3, 5, 10} specified by numerical labels. The dashed 
line indicates the maximum entanglement achievable in the 
asymptotic limit M — ► oo, squeezing constant r — 0.5 and 
Co = 1. 
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FIG. 5: The same as Fig. 4 but the logarithmic negativity En 
of the two-mode atomic state is plotted as a function of the 
resonant optical depth d of the atomic ensemble. The various 
lines correspond to different number of Gaussification steps 
M = 1 (solid line), M = 2 (dashed line), M = 3 (dashed- 
dotted line) and M = 4 (dotted line). Shown are results 
for protocol with variable coupling strength n (a), and for a 
protocol with fixed n — 1 (b). The parameters read r — 0.5, 
n = 0, and Co = 1. 



see that each step of the protocol increases entanglement 
of the two memories. For a fixed number of steps the 
amount of created entanglement monotonically decreases 
with increasing amount of thermal noise n. 

So far we have considered an idealized scenario with 
no decoherence of the atomic memory. We now inves- 
tigate the decoherence effects caused by finite resonant 
optical depth d and nonzero atomic depumping r/. In 
particular, we compare the Gaussification protocol sug- 
gested in this paper with the simple scheme with fixed 
coupling strength as described by Eq. ((6|). These two 
protocols are both used to prepare two-mode squeezed 
vacuum state of two atomic ensembles by mapping many 
copies of a two-mode squeezed vacuum state of light onto 
atoms. In a limit of infinite d and zero 77, both protocols 
are completely equivalent and create the same amount 
of entanglement after M Gaussification steps. However, 
the atomic decoherence influences the two schemes in a 
different way, because they exhibit different K^ ot . We 
have numerically calculated covariance matrix of a two- 
mode atomic state after M steps of each protocol. In 
the numerical simulation, the decoherence map (J8J) with 
appropriate depumping factor rjj = Kj/d is applied to a 
covariance matrix of the atomic state after each interac- 
tion of atoms with light. The results are reported in Fig. 
5 where we plot the logarithmic negativity of the atomic 
state as a function of d for several steps of the protocol, 
M = 1,2,3,4. A comparison of Figs. 5(a) and 5(b) re- 
veals that the protocol exploring switching between two 
different QND couplings and variable coupling strength 
can generate substantially more entanglement than the 
protocol with fixed k. In this latter case already the sec- 
ond step can be in fact detrimental and lead to reduction 
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FIG. 6: The same as Fig. 4 but logarithmic negativity of 
the atomic state En after 2 (a) and 4 (b) steps of the proto- 
col is plotted as a function of homodyne detection efficiency 
?7hd- The results are shown for protocol with variable cou- 
pling strength (solid line) and for a protocol with fixed k = 1 
(dashed line). The parameters read r = 0.5, n — 0, and 
Co = 1. 



of atomic entanglement due to decoherence if d is small 
enough. 

Besides atomic depumping, another effect that can 
negatively influence the Gaussification procedure is im- 
perfect homodyne detection that is used to measure the 
output light quadrature. The quadrature operator de- 
tected by an imperfect homodyne detector with detection 
efficiency ?7hd reads p L , e $ = v^hdPl + V 1 = VhbPl, noise 
where pi, noise represents quadrature of an auxiliary vac- 
uum state. In the feedback, we can compensate for 
ijhd < 1 by re-scaling the gain, g — > g / ^/?7hd ■ The noise 
added to an atomic quadrature by feedback is then given 

by 



1HD 



-PL, 



This noise increases both with de- 



creasing efficiency of homodyne detection and with in- 
creasing feedback gain. The protocol with fixed k speci- 
fied by Eq. (|6]) is less sensitive to this decoherence mech- 
anism, because there the gain scales with number of steps 
M as 1 /M while for the scheme with varying k the scaling 
is I/Af 1 / 2 . This is confirmed by numerical calculations 
whose results are shown in Fig. 6. In order to unam- 
biguously assess the influence of ?7hd, we have neglected 
other decoherence mechanisms and considered the limit 
of d —> oo. We can see that the entanglement generated 
in the atomic memory decreases with decreasing ?7hd- 
However, both protocols are rather robust and atomic 
entanglement can be generated even if ?7hd < 50% while 
typical efficiency of homodyne detection can exceed 90%. 
Note also that in the considered example the difference 
between the two protocols becomes non-negligible only 
for homodyne detection efficiency less than 80%. 



IV. ENTANGLEMENT PURIFICATION 

As a final application of the Gaussification protocol we 
consider purification and distillation of non-Gaussian en- 



tanglement of light beams coupled to atomic memories 
placed at two distant locations A and B. The protocol 
is again based on local interactions of parts of entangled 
states with atomic memory. Additionally, we also make 
use of conditioning on the outcomes of homodyne mea- 
surements on light beams after they interacted with the 
atomic memory [3, E3) ■ The homodyne detection thus 
plays a dual role here: on one hand it provides the in- 
formation required for a feedback onto atoms and on the 
other hand it heralds the success of purification. In par- 
ticular, the scheme succeeds if all measurement outcomes 
Qj are sufficiently close to zero, \Qj\ < Qt, where Qt is 
some threshold. In what follows we will consider limit of 
a very narrow acceptance window which can be mathe- 
matically modeled as conditioning on Qj =0. Note that 
the acceptance condition must be simultaneously satis- 
fied at both locations A and B. The protocol therefore re- 
quires classical communication between A and B which is 
used to establish success or failure. Although the present 
scheme involves quantum memories, it is not a quantum 
repeater because the memories do not improve the scaling 
of the success probability of the protocol. However, the 
quantum memories do enable purification and Gaussifi- 
cation of many copies of entangled states of light emitted 
at different times (l9j |. 

As an explicit example we consider purification of 
phase-diffused two-mode squeezed states which has previ- 
ously been successfully demonstrated experimentally for 
traveling light beams [l7j]. Phase diffusion arises when 
the entanglement is distributed over quantum channels 
with fluctuating optical length. We assume that the 
source is located at a center between A and B such that 
both modes suffer from random phase shifts 4> A and 4>b ■ 
The phase noise converts the initial pure entangled Gaus- 
sian two-mode squeezed state \^tms)ab into a mixed 
non-Gaussian state 



Pab = 



where ^ab 



{U A ®U B )^ AB {u\®ul)P{cj> A )P^ B )d^ A d(j )E 

(30) 



I* 



TMS 



X* 



TMS I 



u, 



exp(—i(f>jiij), 



denotes photon number operator of mode j, and P{4>j) 
is the probability distribution of random phase shift. In 
what follows we shall assume Gaussian distribution of 
random phase shift [l3|, H3] j 



1 



V27 



: exp 



2a 2 



(31) 



where a quantifies the strength of fluctuations. Note 
that Gaussian distribution of cj> gives rise to highly non- 
Gaussian p AB - For strong enough phase noise, the Gaus- 
sian entanglement is completely lost and the two-mode 
state does not exhibit any quadrature squeezing. In such 
case the deterministic Gaussification protocol discussed 
in previous section is useless, because it would produce a 
separable state of the two memories at stations A and B. 
In contrast, the protocol augmented by conditioning on 
outcomes of homodyne detection is capable of converting 
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the initially non-Gaussian entanglement into Gaussian 
one. We quantify the performance of the protocol by 
logarithmic negativity En of the distilled state p, purity 
of the distilled state V = Tr(p 2 ), and total variance 

1 = i (((Ax A - Ax B ) 2 ) + ((A PA + A PB f)) . (32) 

The state is entangled if X < 1 . We also calculate Gaus- 
sianity of the state, G, which is defined as fidelity of the 
state with a Gaussian state with the same covariance ma- 
trix and displacement. It holds that G < 1 and G = 1 
only if the state is Gaussian. 

The results of numerical simulations are shown in Fig. 
7, which shows the total variance, logarithmic negativity, 
purity and Gaussianity of the state in atomic memories 
after two and twenty steps of the protocol. For compari- 
son, the thick solid gray lines indicate values of the con- 
sidered quantities attained by the de-phased light state. 
Each successful step of the protocol increases purity and 
Gaussianity and after 20 steps of the protocol the state is 
very pure and very close to Gaussian state for all a < 1 . 
Each successful measurement step also increases logarith- 
mic negativity of the atomic state and decreases its total 
variance. For large a the de-phased state of light exhibits 
X > 1, i.e. no Gaussian entanglement. The purification 
reduces the total variance such that X < 1, hence it con- 
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Psucc Qt 

FIG. 8: Entanglement purification of phase-diffused two- 
mode squeezed states with non-zero acceptance threshold 
Qt > 0. Figure (a) shows the total variance of atomic state 
as a function of the overall success probability of the protocol. 
In figure (b) we plot the dependence of overall success prob- 
ability on acceptance threshold Qt- The results are shown 
for a two-step protocol (solid line) and for a four-step proto- 
col (dashed line). The parameters read a = 1, r = 0.5 and 
Co = a/Vtt. 
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FIG. 7: Gaussification of phase-diffused two-mode squeezed 
vacuum into two distant atomic memories. The figure shows 
total variance (a), logarithmic negativity (b), purity (c), and 
Gaussianity (d) of the atomic state after 2 (dashed line) and 
20 (dotted line) steps of the protocol as functions of phase fluc- 
tuations strength a. For comparison, solid gray line shows the 
values corresponding to the initial phase-diffused two-mode 
squeezed state of light. The atomic ensembles are initially in 
vacuum state, the two-mode squeezing constant r — 0.5 and 
Go = J2j^. 



verts non-Gaussian entanglement into Gaussian one. In 
the present example, the entanglement stored in atomic 
memories does not exceed the initial entanglement of the 
de-phased light beam c.f. Fig. 7(b). One reason for 
this saturation effect is that the two atomic ensembles 
are initially in a separable state, in contrast to entan- 
glement distillation protocols with light beams where all 
copies are initially in the same entangled state [l3T[l5l |. 
Secondly, the procedure somewhat differs from the pre- 
viously studied Gaussification-based entanglement distil- 
lation schemes due to switching between measurements 
of amplitude and phase quadratures. 

An important parameter of the entanglement purifi- 
cation protocol is its success probability P SUC c- In order 
to assess how the performance of the scheme scales with 
^succ we have performed numerical simulations of a full 
protocol with non-zero acceptance threshold Qt and de- 
termined the dependence of total variance X of the state 
of two atomic ensembles on the overall success probabil- 
ity of the protocol. The calculations were performed for 
two-step and four-step scheme. Results reported in Fig. 
8 clearly demonstrate that the probabilistic entanglement 
purification protocol can significantly outperform the de- 
terministic scheme for an overall success probability of 
10% or even higher. This is in full agreement with pre- 
vious theoretical and experimental results on distillation 
of phase diffused squeezed and entangled states of light 
where als o g ood performance for large success rates was 
observed [13, El]. 
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V. CONCLUSIONS 

In summary, we have proposed and analyzed a protocol 
for Gaussification of states of M traveling light beams in 
an atomic memory. The protocol is designed for a mem- 
ory that interacts with light via quantum non-demolition 
coupling. The scheme involves repeated switching be- 
tween xaPl and pa%l couplings, homodyne detection 
of output light and feedback on atoms. In contrast to 
other proposals, the total coupling strength required by 
this procedure scales only logarithmically with the num- 
ber of Gaussified modes and no presqueezing operation 
is required on atoms or light. All components of the 
proposed protocol have already been successfully exper- 
imentally demonstrated in the past [3, H| so small-scale 
demonstrations with M = 2 or M = 3 should be feasible 
with present-day technology. If combined with condi- 



tioning on the outcomes of homodyne detection on out- 
put light, the present protocol can be used to purify 
and Gaussify non-Gaussian mixed entangled states. This 
procedure can thus for instance increase performance of 
certain entanglement-based continuous variable quantum 
key distribution schemes that require high purity and 
Gaussianity of shared entangled state for optimal perfor- 
mance. 
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